Spin dephasing and pumping in graphene due to random spin-orbit interaction 
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We consider spin effects related to the random spin-orbit interaction in graphene. Such a random 
interaction can result from the presence of ripples and/or other inhomogeneities at the graphene 
surface. We show that the random spin-orbit interaction generally reduces the spin dephasing 
(relaxation) time, even if the interaction vanishes on average. Moreover, the random spin-orbit 
coupling also allows for spin manipulation with an external electric field. Due to the spin-fiip 
interband as well as intraband optical transitions, the spin density can be effectively generated by 
periodic electric field in a relatively broad range of frequencies. 

PACS numbers: 72.25.Hg,72.25.Rb,81.05.uc,85.75.-d 



I. INTRODUCTION 

Graphene is currently attracting much attention as a 
new excellent material for modern electronics^"— . The 
natural two-dimensionality of graphene matches per- 
fectly to the dominating planar technology of other semi- 
conducting materials, and correspondingly gives the way 
to creating new hybrid systems. However, the most strik- 
ing properties of graphene are not directly related to its 
two-dimensionality. Due to the bandstructure effects, 
electrons in pure graphene can be described by the rela- 
tivistic Dirac Hamiltonian, leading to the linear electron 
energy spectrum near the Dirac points. As a result, the 
electronic and transport properties of graphene are sig- 
nificantly different from those of any other metallic or 
semiconducting material;^"— except (to some extent) its 
parent material - the clean graphitei^ 

It has been also suggested that graphene may have 
good perspectives as a new material for applications 
in spintronicsi^"— The intrinsic spin-orbit interaction 
in graphene is usually very small, and therefore one 
can expect extremely long spin dephasing (relaxation) 
timeJ^"— Thus, spin injected to graphene, for instance 
from ferromagnetic contacts, can maintain its coherence 
for a relatively long time. Experiments demonstrate 
spin relaxation times for various graphene-based systems 
spanned over several orders of magnitudoi^r— with some 
of them being much shorter than expccted i^^i^^ The rea- 
son of this contradiction is not quite clear, and several 
different explanations of these observations have been 
already put forward i^"— In this paper we present an- 
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other model based on the random Rashba spin-orbit 
interaction.— Physical origin of such random spin- 
orbit interaction can be related to the ripples existing 
at the surface of graphene^Zr— and/or to some impuri- 
ties adsorbed at the surface, which randomly enhance the 
magnitude of spin-orbit coupling as compared to that in 
the clean graphene<^ 

One of the key issues in spintronics (including 
graphene-based spintronics) is the possibility of spin ma- 
nipulation with an external electric/optical field. This 
includes spin generation, spin rotation, spin switching, 
etc. Here we consider the possibility of spin pumping in 
graphene using the idea of combined resonance in systems 
with Rashba spin-orbit interactioni^i^ The possibility 
of spin manipulation using optical excitation^ is based 
on various mechanisms of spin-orbit interaction in semi- 
conductor systems. In particular, the spin polarization 
appears in systems with regular spin-orbit coupling, sub- 
ject to periodic electric fieldi^"— It has been shown re- 
cently, that the random spin-orbit interaction also can be 
applied to generate spin polarization in symmetric semi- 
conductor quantum wcUsi^ In this paper we show that 
similar method can be used to generate spin polarization 
in graphene with random Rashba spin-orbit interaction. 
To do this we analyze the intensity of optically-induced 
spin-flip transitions assuming two-dimensional massless 
Dirac model of electron energy spectrum in graphene, 
and calculate the magnitude of spin-polarization induced 
by the optical pumping. 

The paper is structured as follows. In section 2 we 
describe the model Hamiltonian assumed for graphene. 
Spin dephasing due to the random spin-orbit Rashba cou- 
pling is calculated in section 3. In turn, spin pumping by 
an external electric field is considered in section 4. Final 
conclusions arc presented in section 5. 
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II. MODEL 

To describe electrons and holes in the vicinity of Dirac 
points we use the model Hamiltonian "Hp which is suffi- 
cient when considering the effects related to low-energy 
electron and hole excitations. We also include the spin- 
dependent perturbation in the form of a spatially fluc- 
tuating Rashba spin-orbit interaction, 'Hso- Thus, the 
system Hamiltonian can be written as (we use system of 
units with h = 1) 



T~L — Ho + Hsoi 
Ho = VT ■'k, 

^° V \ dkx ^ dky ^ 
= A(r) {r^ay - Tya^;) , 



(1) 
(2) 



(3) 



where v is the electron velocity, A(r) is the random spin- 
orbit parameter, r = (x, y) is the two-dimensional coor- 
dinate, and T and cr are the Pauli matrices acting in the 
sublattice and spin spaces, respectively.— Equations ([2]) 
and ([3]) show that spin-orbit coupling can be described 
by a conventional Rashba Hamiltonian, proportional to 
VxOy — VyUxi where the velocity components Vx, Vy are, in 
general, obtained with the unperturbed Hamiltonian "Ho- 
lt is well known that there is an intrinsic (internal) spin- 
orbit coupling in graphene, which is related to relativistic 
corrections to the crystal field of the corresponding lat- 
tice. In addition, a spatially uniform Rashba field can 
be induced by the substrate on which the graphene sheet 
is located. The reported results suggest that these inter- 
actions, which can be considered as independent sources 
of spin relaxation, are cither very weaki^i^ or do not 
contribute to the dcphasing rate by symmetry reasonsi^ 
Therefore, we will neglect them in our considerations and 
will briefly discuss their role in the following. 

The Schrodinger equation, (Ho — e) ?/'k = 0, for the 
pseudospinor components of the wavefunction ip^ is 



Xk 



0, 



(4) 



kx ± iky. The normalized solutions corre- 



where k± 

sponding to the eigenstates Sk = ±ufc of Hamiltonian Ho 
can be written in the form 



V'kfT±(r) 



V2 



la)±-i|2a) 



(5) 



where the ± signs correspond to the states in upper and 
lower branches, respectively. 

We assume that the average value of spin-orbit inter- 
action vanishes, while the spatial fluctuation of A(r) can 
be described by the correlation function F {r — r') of a 
certain form, 

(A(r))=0, (6) 
C(r - r') = (A(r) A(r')) = (A^) Fir- r') . (7) 



When calculating spin dephasing, one can consider 
only the electron states corresponding to the upper 
branch (conduction band) of the energy spectrum, ek ~ 
vk. The spin-flip scattering from the random potential 
determines the spin relaxation in this particular band, 
while from symmetry of the system follows that spin de- 
phasing in the lower (valence) band is the same. The 
intraband matrix elements of the random spin-orbit in- 
teraction ([3]) in the basis of wavefunctions (O for the 
conduction band form the following matrix in the spin 
subspace 



Vkk' = A 



kk' 





ik'_^/k' 



~ik^/k 




(8) 



where Akk' is the Fourier component of the random spin- 
orbit coupling. Since scattering from the random spin- 
orbit potential is elastic, only the intraband transitions 
contribute to the spin relaxation. 



III. SPIN DEPHASING 

To demonstrate how the random spin-orbit coupling 
works in graphene and how its effects can be observed 
in experiment, as well as to compare graphene and con- 
ventional two-dimensional semiconductor structures, we 
calculate in this Section the corresponding spin dephas- 
ing time. For this purpose we use the kinetic equation for 
the density matrix (Wigner distribution function) j^i^ 



9pk 

dt 



Stpk- 



(9) 



The collision integral St pk on the right-hand side of this 
equation is due to the spin-flip scattering from the ran- 
dom spin-orbit interaction, 

Stpk = (2Vkk'Pk'Vk'k - Vkk'Vk'k/Ok - PkVkk'Vk'k) 

k' 

x<5(ek-ek')- (10) 
We assume the following form of the density matrix; 



Pk = Pak + SkO-z 



(11) 



where the flrst term pofc corresponds to the spin- 
unpolarized equilibrium state. On substituting ([5]) and 
pT]) into Eq. ^ we flnd 



Stpk = -27r(72^C(q)(sk + Sk') <^(ek-ek')' 



(12) 



where C(q) = (A|j.,) and q = k' — k. Assuming that 
Sk does not depend on the point at the Fermi surface we 
obtain 



St pk = - 



Sk 



2k 



^Ak?- — g' 



■.dq. 



(13) 
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FIG. 1: Spin dephasing time as a function of the characteristic 
range R of the random spin-orbit fluctuations. 



For dcfiniteness, we assume that the characteristic spa- 
tial range of the random spin-orbit fluctuations is R, and 
take C(q) in the following form: 



C(q) = 27r(A2)i?2 



-qR 



satisfying the normalization condition 



C(q) 



(2^ 



(14) 



(15) 



The resulting spin relaxation rate is not strongly sensi- 
tive to the shape of the correlator. However, the ap- 
plicability of the approach based on Eg. p^ depends on 
the ratio of electron mean free path £ to i?, being valid 
only ii £/R ^ 1. In such a case, typically realized in 
graphene, the electron spin experiences indeed random, 
weakly correlated in time fluctuations of the spin-orbit 
coupling. In addition, we can safely neglect the effect of 
the random spin-orbit coupling on the momentum relax- 
ation rate. Finally, for the spin dephasing time we obtain 
the following expression: 



ink 



2kR 



(A2) i?2 [U2kR) - Lo(2kR)] . 

V 



(16) 



where /o(a;) and Lo{x) are the modified Bessel and Struve 
functions of zeroth order, respectively. 

In the limiting semiclassical {kR 3> 1) and quantum 
{kR ^ 1) cases we find 



1, kR > 1, 



irkR, fci? < 1. 



(17) 



For kR ^ 1, the result in Eq. p7|) can be interpreted as 
the special realization of the Dyakonov-Perel' spin relax- 
ation mechanism. To see this, we note that the electron 

1 /2 

spin rotates at the rate $7 ~ (A^) , with the preces- 
sion direction changing randomly at the timescale of the 
order of time that electron needs to pass through one do- 
main of the size R, i.e., tr R/v. The resulting spin 



relaxation rate l/r| is of the order of H.'^tr. It is worth 
mentioning that at given spatial and energy scale of the 
fluctuating spin-orbit field, the decrease in the electron 
free path and in the momentum relaxation time leads to 
the decrease in the spin dephasing rate: if £ <C i?, elec- 
tron spin interacts with the local rather than with the 
rapidly changing random field, and spin relaxation rate 
becomes of the order of (A^) r, where r is the momen- 
tum relaxation time. This agrees qualitatively with the 
observations of Ref.jl^, however, a quantitative compar- 
ison needs a more detailed analysis. 

Taking examples with typical values v = 10^ cm/s, 
R ^ 50 nm, and (A^) ~ 500 /xeV^, similar to what can 
be expected from Ref.[l2l, we obtain t| less than or of 
the order of 10 ns. As one can see from Ea. p?|) . the 
spin relaxation for small kR <^ 1 is suppressed, as can 
be understood in terms of the averaging of the random 
field over a large 1/k^ ^ R^ area. The full fc-dependence 
in Ea. (jl7p implies that the spin dephasing rate is pro- 
portional to n^/^ at low carrier concentrations n and 
is independent of n at higher ones. Therefore, at the 
charge neutrality point, where n = 0, the spin relaxation 



vanishes, in agreement with the observations of Ref. [21 



Moreover, our approach qualitatively agrees with the in- 
crease in the spin relaxation time in the bilayer graphene 
compared to the single layer one^i: the transverse rigidity 
of the bilayer can be larger, thus suppressing formation 
of the long-range ripples, and, as a result, the random 
spin-orbit coupling. 

Equation pT)) shows that, as far as the spin dephasing 
is considered, the only difference between graphene and 
conventional semiconductor a^^'^^i'^^ is related to the fact 
that the electron velocity is constant for the former case 
and is proportional to the momentum in the latter one. 
As we will sec in the next Section, this difference becomes 
crucial for the spin pumping processes. 

The dependence of spin relaxation time, calculated 
from Eq. (jl6|) as a function of the characteristic domain 
size R of the random spin-orbit interaction is presented 
in Fig. 1, were Tsq is defined as t^q^ = 8(A^) /vk. The 
curves corresponding to different values of /c > 10^ cm~^ 
would be practically indistinguishable in this figure. 

Here several comments on the numerical values of spin 
relaxation are in order. The values observed in experi- 
ments on spin injection from ferromagnetic contact o^^i^^ 
are of the order of 10"^*^ s, two orders of magnitude less 
than our estimate which does not take into account ex- 
plicitly the role of the Si-based substrates. The effect of 
the Si02 substrate, including the contributions from im- 
purities and electron-phonon coupling, was thoroughly 
analyzed in Ref. [23| . However, the obtained dephas- 
ing rates were well below the experimental values and 
also below the estimate obtained here, leading the au- 
thors of Ref.jl^ to the suggestion of an important role 
of heavy adatoms in the spin-orbit coupling.— On the 
other hand, it was shown that the spin dephasing rate 
can be strongly infiuenced at relatively high temperatures 
by the electron-electron collisions However, including 
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these collisions does not bring theoretical values closer to 
the experimental ones. 

The discrepancies between theory and experiment and 
between experimental data obtained on different systems 
call for a more detailed analysis of the experimental sit- 
uation, including the dependence of spin relaxation time 
on the device functional properties and the experimental 
techniques applied. 

IV. SPIN PUMPING 

Let us consider now spin pumping by an external elec- 
tromagnetic periodic field corresponding to the vector 
potential A{t). Wc assume that the system described by 
Eqs. (l)-(3) is additionally in a constant magnetic field 
B. For simplicity, we consider the Voigt geometry with 
the field in the graphene plane, so that the effects of Lan- 
dau quantization are absent. Thus, the Hamiltonian "Ho 
includes now the constant field and can be written as 

no = vT-k + Aa^, (18) 

where 2 A — gfisB is the spin splitting and the magnetic 
field is oriented along the a;— axis (the electron Lande fac- 
tor for graphene is g = 2). The induced spin polarization 
is opposite to the direction of the magnetic field. 

Since we are interested in real transitions in which the 
energy is conserved, we consider interaction with a single 
component of the periodic electromagnetic field, A(t) = 
Ae~'"*, which enters in the gauge-invariant form 

nA = — ^-A{t)=--T-A{t), (19) 
c ok c 

and in the following we treat the term Ha as a small 
perturbation. 

The absorption in a periodic field (probability of field- 
induced transitions in unit time) can be written as 

/de 
— nAGi,is + Lj)nAGk{e), (20) 

K ^ 

where Gk(e) is the Green function. In the absence of 
spin-orbit interaction, the absorption ()20p does not in- 
clude any spin-fiip transitions. We can account for the 
spin-orbit interaction ([3|) in the second order perturba- 
tion theory, including the corresponding matrix elements 
as shown in Fig. 2. This means that wc do not con- 
sider perturbation terms as the self-energy within a single 
Green function assuming that they are already included 
in the electron relaxation rate. 

Thus, in the second order perturbation theory with re- 
spect to the random spin-orbit interaction Jiso we obtain 

I{uj) = Rc Tr ^ / 1^ Ha Gi{e + u) H^^' [e + u) 

kk' ^ 

>^nAGl,{e)H^i'^Gl{e)^ (21) 




FIG. 2: (Color online.) Feynman diagrams for the light ab- 
sorption in the second order perturbation theory with respect 
to the random spin-orbit interaction. The coupling and elec- 
tromagnetic vertices are shown as white and filled circles, re- 
spectively. 




FIG. 3: (Color online.) The energy spectrum and indirect 
spin-flip transitions in graphene in a uniform magnetic field 
B. Long and short arrows correspond to the interband and 
intraband transitions, respectively. 

where Green's function G^(£) = diag |G^.j.(e), G|^j^(£)| 
corresponds to the Hamiltonian "Ho of Eq. (fT8|) . 

g,o / ^ £ + i;T-k + Acr + /^ 

[e - eika +1-1 + iS sgn e) (e - e2ka + i-i + iS sgn e) ' 

(22) 

with (T = ±1 corresponding to the spins oriented along 
and opposite to the a;— axis, respectively, £(i,2)fe(T ~ 
±v k + Act, S being the half of the momentum relaxation 
rate, S = 1/2t, and n denoting the chemical potential. 

Diagrams in Fig. 2 show the qualitative difference 
between the graphene and semiconductor quantum well 
with respect to the effects of random spin-orbit coupling. 
In semiconductors, the external electric field is explicitly 
coupled to the anomalous spin-dependent term in the 
electron velocity, which is random, and therefore the di- 
agrams describing the corresponding transitions include 
only two Green functions. In graphene, due to the ab- 
sence of randomness-originated term in the Hamiltonian 
Ha, four Green functions are required to take into ac- 
count the random contribution of the spin-orbit coupling. 
This situation, in some sense, is more close to what is ob- 
served in the conventional kinetic theory of normal metal 
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conductivity, where the couphng to the external field 
does not depend on the randomness explicitly, and the 
additional disorder efi^ects appear due to the self-energy 
and/or due to the vertex corrections, as in Fig. 2. 

Upon calculating contributions from the diagrams of 
Fig. 2, one finds the total rate of spin-fiip and spin- 
conserving transitions due to the random spin-orbit cou- 
pling in the form 



Jrso(w) = 



ReTr5]5^.2C(q) / ^(r-n^) 

crcr' kk' •' ^ 

xGl,,,{e)T+Gi,{e), (23) 



where t± = ± iTy, and describes the direction of 
A. In the following we assume linear polarization of 
light, A = {A, 0). After calculating the trace and in- 
tegrating over energy e in Eq. one finds a rather 
cumbersome expression (see Appendix) consisting of sev- 
eral terms, each of them corresponding to transitions be- 
tween certain branches of the spectrum (Fig. 3). For 
definiteness, we locate the chemical potential in the va- 
lence band. Correspondingly, only the transitions from 
the bands (2 t) and (2 |) to the unoccupied states in 
the bands (1 f), (1 i), (2 t) and (2 |) are possible. We 
will concentrate on the optically induced spin-flip transi- 
tions contributing to the optically-generated spin pump- 
ing. Hence, we do not consider spin-conserved transitions 
contributing to the usual Drude conductivity!^"— 



A. Interband spin-flip transitions 

Let us consider first the spin-fiip transitions from the 
valence to conduction bands, such as k2 — >■ kj^. For con- 
venience we introduce the parameter kj = (k, cr^)^ , de- 
scribing the momentum and spin projection for an elec- 
tron in the subband j. The corresponding expression for 
the transition rate can be obtained from the equations 
presented in the Appendix, and considerably simplified 
by taking into account that: (i) in the nonsingular terms 
~^k, can be substituted by w, (ii) in the semiclassical 
limit kR 3> 1 the energy change due to the change in the 
momentum is small compared to 1 /r, and (iii) the photon 
energy is much larger than the characteristic low-energy 
scale parameters, i.e., oj » 1/r, and w » A. We mention 
that the linear in A terms have to be kept despite A <C a; 
since the resulting spin pumping rate, determined by the 
contributions of both initial spin states, is linear in A. 
The expressions for the parameters K2a-{'-^), which deter- 
mine the spin-fiip rate as introduced in the Appendix, 
Ea. (|A2p . can be then simplified, and as a result one ob- 
tains the following formula for the spin-flip rate in the 
relevant frequency domain: 



ha^i.'iu^) = 4na^- ^ G (q) Ifie^J - fie^, ] 
c OJ ^ — ' L ^ 1 

kk' 
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FIG. 4: (Color online.) The parameter K2i{i^) for the rate 
of the spin-flip interband transitions in the high-frequency 
domain (a) and K2f (b) for difi^erent values of the range pa- 
rameter R describing characteristic size of the fluctuations in 
spin-orbit interaction and (A^) = 100 /xeV^ . 



fc^ — kq cos ip 



UJ — Er' 



- UJ — £-7-1 



k2 k. 



.(24) 



The transitions are constrained due to the (5-function 
in Eq. (j24p corresponding to the energy conservation 
with the change in the momentum q (here we use (p = 
cos-i(k, q)) 



\lj — vk — 2crA| 5{lp — ipo) 
v'^kq I siniy9o| 



(25) 



where ipo is a solution of the equation 



vk + 2aA-UJ + v^/k^ - 2kqcoaip + q^ = 0, (26) 

which gives us the condition for a minimum value of mo- 
mentum in Eq. (j26p . vkmin = w — 2ctA. The energy con- 
servation determines the angle ipo between the vectors k 
and q as 



cos^o = | 



uj - 2a A ( ^ uj- 2a A 



vq 



2vk 



(27) 



The usual condition of | cos (po\ < 1 leads to the following 
restrictions in the integration over q in Eq. (|24p : 

a) if uj/2- aA < vk < uj~2aA then vk - \vk^2aA- 
Lo\ < vq < vk + \vk + 2a A — w|, 

b) if a;/2-CTA > vk then -vk^\vk^2aA-ijj\ < vq < 
vk + \vk + 2a A - uj\. 

Accounting for all these conditions, the integral over 
k and q in Eq. ^2^ can be calculated numerically. We 
used the following parameters: /i = —3 meV, correspond- 
ing to the Fermi momentum fcj? « 5 x 10^ cm~^ and 
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FIG. 5: (Color online.) The parameter K2{(jj) = Kiiiuj) — 
K2^{uj) for the spin injection rate for different R. Other pa- 
rameters are the same as Fig.4. 



hole concentration w 4 x 10® cm"^, A = 1 meV, and 
1/r = 1 meV. These parameters correspond to a rather 
clean graphene and strong Zeeman splitting of the bands 
(magnetic field of 17 T). For numerical accuracy, the cal- 
culations were performed with exact Eas. (jA4p and (jA6p . 
The results of numerical calculations for the interband 
spin-flip transitions are presented in Fig. 4. Figure 5 
presents the quantity describing the total spin injection 
rate for the interband transitions in this frequency do- 
main, K2{uj) = K2i{uj) — K2-\-{uj). The positive sign of 
K2{uj) corresponds to the fact that the absolute value of 
spin density decreases due to the pumping. The peaks 
correspond to the absorption edge of direct optical transi- 
tions, where the transition probability rapidly increases, 
and the increase at large frequencies is due to the linear 
energy dependence of the density of states. 



Intraband transitions in the hole subband 
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FIG. 6: (Color online.) The parameter K2i{i^) for the rate 
of the spin-flip intraband transitions in the low-frequency do- 
main (a) and K^tiuj) (b) for different values of the range pa- 
rameter R. Other parameters are the same as Fig.4. 
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Using the general formulas (jA4p and (|A6[) we can 
write down the expression for the spin-flip intraband 
transitions within the valence band in a low-frequency 
region uj <^ \fi\. Such transitions are associated with 
a relatively large change of the momentum in the ab- 
sorption process, and one can expect that they give a 
smaller transition rate. Upon taking into account that 

^ we find 



-ki 



-2/1 ^ ui, and 



-in a- 



E 

kk' 



^^k2^k' ^^'^ ^ kqcosip) 



6 (eT-^ +0; - a 



(e^ +w-e^O(e^ -e^O 



,(28) 



The (5-function can be presented in the form of Eq. (|25p , 
where (po is a solution of the equation 



vk + 2a A 



■ UJ — V 



y/k^ -2kqcosip + q^ = 0. (29) 



FIG. 7: (Color online.) The parameter 7^2 (cj) = K2i{to) — 
K2f{i^) for the spin injection rate for different R in the low- 
frequency domain. Other parameters are the same as Fig.4. 
Unmarked line corresponds to i? = 30 nm. 



Equations and (pS)) are similar, with 

the important difference in the sign in front of 
\J k"^ — 2kq cos tf + . In the case of interband tran- 
sitions this corresponds to the transitions occurring 
at w sa 2vk, while the intraband transitions occur at 
lower frequencies determined by the possible momentum 
transfer due to the randomness of the spin-orbit cou- 
pling. The solution exists only for vk > u) ~ 2crA and 
gives Ea. (|27|) . However, the condition |cos(^o| < 1 leads 
here to a different restriction. Since only the condition 
Lo — 2(tA < 2vk is consistent with vk > co — 2aA, the 
momentum q should be then in the single range of 
vk — \vk + 2a A ~ uj\ < vq < vk + \vk + 2a A — w|, in 
contrast to the case of interband transitions. 

The results of calculations for the intraband transitions 
arc presented in Fig. 6. The intensity of such processes 
is relatively small compared to the interband transitions. 
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and they can be seen only at low photon energies. Fig- 
ure 7 corresponds to spin injection rate by the intraband 
transitions, K2{u)). 

V. CONCLUSIONS 

We have considered certain spin effects associated with 
random spin-orbit interaction in graphene. First, we have 
calculated the corresponding spin relaxation time, and 
believe that this mechanism can be dominating when the 
amplitude of fluctuations in spin-orbit interaction is large 
enough. This may happen in the presence of surface rip- 
ples with short wavelengths. The other possibility can be 
related to the absorbed impurities at both surfaces of a 
free-standing graphene. One can expect especially strong 
random spin-orbit coupling for heavy impurity atoms. 

The second effect concerns the possibility of spin 
pumping by an external electromagnetic field. The re- 
sults of our calculations show that graphene can be used 
as a material, in which the electron spin density can be 
generated by the optical pumping. The mechanism of 
pumping here is related to the spin-flip transitions asso- 
ciated with the random Rashba spin-orbit interaction. 
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Appendix A: Formula for the absorption rate 

Using ([^ and after calculating the trace and 

integrating over energy e, we find out that spin- flip 
processes can be characterized by the initial state as 
/il(a;), J24,(a;), and /2t('^) with the corresponding 
transition rate /j[t(w) defined as 



hp A 

I.^ic) ^ ^^K.^ico), (Al) 

f \ T /" 4^/ Jj„{(i,uj) (fk(Pk' 

K,A.) ^ -ImJ V dd) ^^^^ , (A2) 



describing transitions from the band and spin states cor- 
responding to the subscript ja. Here Jj^{q,Lu) have the 
form: 



, , , j.r , „ , M (slki- +£ikiuj){e ^kq COS if) 
Ji;(q,w) = Ifieiki) - /(eifei + 



1 



{eiki +UJ - eik'^ + iS) {eiki + w - e2ki + i5) (eiki + - e-ik't + «<^) 



'/2i(q,w) [fie2ki.) - f{e2ki+^)] ' ^''^ 



{£2ki — £2fe't)(£2fc4. — £lki){£2kl — ^Ifc't) 
1 



{£2ki +0J - e2k't + iS) {£2ki + w - siki + iS) {e2ki + - sik't + i5) 
Jit(q,w) = [/(eu-t) - /(eifef + ^^)] 



(e^j.^ - + eik\ijj){k'^ - kqcosip) 



(^ifct ^ £ifc'i)(£ifet ~ e2fet)(^ifet ~ ^2k'\) 
1 



(eifc-t- + w - eifc'^ + iS) (eifct + ^ - £2fet + i5) (ei^t + - £2k'i + iS) 
(sjfet - + £2k^oj){k'^ - kq cos (f) 



J2t(q,w) ^ [/(e2fet) - /(£2fet + 



(£2fet ~ £2fe'i)(£2fct ~ eifct)(£2fet ~ £lk'i) 
1 



{e2kt + w - e2k'i + iS) {e2kt + w - eik-t + iS) {s2kt + - sik'i + iS) 

I 



(A3) 



(A4) 



(A5) 



(A6) 



Although the expressions seem to be long, all Jja-{(l,uj) tain energy-difference denominators corresponding to the 
terms have the same simple structure. They con- 
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transitions from initial ja states to all allowed final 
states, and the corresponding resonant terms. Here the 
single allowed spin-conserving transition arising due to 
the Ha term in Ea. (|19p is the momentum-conserving as 
well, while the two transitions caused by random spin- 



orbit term Ti^J^ are off-diagonal both in the spin and 
momentum subspaces, as illustrated in Fig. 2. Taking the 
imaginary part in each of these terms we get 5— functions 
corresponding to the energy conservation for the transi- 
tions to different bands. 
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